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Abstract 

The  existence  of  even  or  odd  solutions  of  Mathieu's  equation  vzith 
period  n  or  2n  depends  on  the  vanishing  of  certain  infinite  determinants.  It 
is  shown  that  these  determinants  can  be  expressed  in  terms  of  the  values  of 
the  even  or  odd  solutions  at  x  =  n/2,  where  x  is  the  independent  variable  . 
The  corresponding  result  for  Hill's  determinant  (which  can  be  expressed  in 
terms  of  the  values  of  certain  solutions  at  x  «  n)  has  long  since  been  found. 
The  method  applied  in  the  present  report  uses  theorems  abouL  the  order  of 
growth  of  the  solutions  both  with  respect  to  x  and  the  parameters.  A  repre- 
sentation of  the  solutions  as  a  finite  Fourier  integral  can  also  be  derived  by 
this  method.  We  indicate  how  far  the  results  can  be  extended  to  the  case  of 
more  general  equations  of  Hill's  type.  The  periodic  solutions  of  the  ijihomo- 
geneous  Mathieu  equations  are  constructed  by  using  an  inversion  formula  for 
certain  infinite  matrices.  The  connection  between  Hill's  equation  and  the 
construction  of  transparent  layers  is  reviewed. 
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Introduction  VJe  shall  write  Mathieu's  equation  in  the  form 

(1)  y"  +  li(co^  +  2t  cos  2x)y  -  0, 

2    2 
where  y"  denotes  d  y/dx  .  The  most  widely  treated  problem  in  the  theory  of 

this  equation  is  the  following  one;  Let  t  be  real.  Find  the  values  of  co  for 

which  there  exists  a  solution  y(x)  of  (l)  such  that  either 

(2)  y(x  +  n)  =  y(x) 
or 

(3)  y(x  +  27t)  =  7(x),  (y(x+  n)  f   y(x)  )  . 

One  nethod  of  finding  a  solution  of  (l)  satisfying  (2)  or  (3)  is  to  satisfy  (l) 
by  an  infinite  series 

oo 
(U)  y  =  H'    ^j.  e^  (i2nx)  , 

n=-oo 

or  by  a  corresponding  series  if  (3)  is  to  be  satisfied.  By  substituting  (h.) 
into  (l)  we  find  a  homogeneous  system  of  infinitely  many  linear  equations  for 
the  c  .  The  determinant  of  this  system  -  aftor  an  appropriate  normalization 
of  the  coefficients  in  the  main  diagonal  -  is  Hill's  determinant,  which  is  de- 
fined explicitly  in  Section  II  and  is  denoted  by  ^  (a),t).  The  vanishing  of 

sin^'TTco  A  is  a  necessary  and  sufficient  condition  for  the  existence  of  a  solu- 
tion of  period  n. 

Another  method  of  dealing  with  the  problem  of  periodic  solution  is 
based  on  Floquet's  theorem.  Using  the  fact  that  the  coefficient  of  y  in  (l) 
is  an  even  continuous  function  of  x  with  period  n,  one  arrives  at  the  follow- 
ing result: 

Let  y, (x)  be  the  solution  of  (l)  defined  by  y, (O)  =  1,  y, (O)  =  0  and  let 
I  .1 

yp(x)  be  the  solution  defined  by  y^(0)  =  0,  yp(0)  =  1.  Then  (l)  has  a  solution 

satisfying  (2)  if  and  only  if 

(5)  y^CTi)  =  1 

and  a  solution  satisfying  (3)  if  and  only  if 

(6)  y^^dO  =  -1  . 
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The  connection  with  the  methods  based  on  Hill's  determinant  is  estab- 
lished by  the  well-known  relation 

(7)  y-^di)  =  1  -  2  Aq  sin^Tico  . 

Recently  SchaefkeE'Ojhas  stressed  the  point  that  for  any  differential 
equation 

(8)  y"  +  q(x)y  -  0, 

in  which  q(x)  is  continuous  for  all  x  and 

(9)  q(x)  =  q(-x),   q(x+7t)  =  q(x)  , 

the  following  elementary  result  holds:  Let  y, (x)  and  yp(x)  be  the  normalized 
solutions  of  (8)  which  are  defined  as  above.  Then  (8)  has 


(10) 


an  even  solution  of  period  n       if  and  only  if  y-|(p')  =  0 
an  odd  solution  of  period  n  "         y^C^)  ■  6 

an  even  solution  of  period  2n      "         y  (^)  =  0 

an  odd  solution  of  period'  2n       "         yp(^)  =  0  . 

The  statements  in  (10)  express  the  periodic  boundary  conditions  in  terms  of 
an  ordinary  boundary  value  problem  for  the  interval  (0,  n/2).  Also,  they 
show  that  the  conditions  expressed  by  (2)  or  (3)  can  be  split  into  two  inde- 
pendent conditions.  This  becomes  obvious  from  the  formulas 

(^^)  yi(Tt)  - 1  =  y[q)  j^q) 

y^(Ti)  *  1  =.  y^(|)  y^q) 

(which,  again,  are  of  an  elementary  nature). 

In  the  case  of  Mathieu's  equation  (l),  the  investigation  of  the  four 
possibilities  listed  under  (10)  is  a  standard  feature  of  the  theory.  Frequently, 
the  problem  of  finding,  for  example,  an  even  solution  of  period  n,  has  been 
treated  by  using  infinite  determinants .  If  such  a  solution  exists  it  can  be 
expanded  in  a  series 

00 

(12)  y(x)  =  ^  g^  cos  2r  x 

r=0 
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On  substituting  (12)  into  (l)  we  find  again  that  a  certain  infinite  determinant 
must  vanish.  We  write  the  resulting  condition  in  the  form 
(13)  -2m  sinnco  C^  =  0  , 

where  C  is  an  infinite  determinant  which  is  given  explicitly  in  equation  (20) 
of  Section  II.  Corresponding  to  the  three  other  oases  listed  in  (10),  three 
more  infinite  determinants  are  defined  in  Section  II;  they  are  denoted  by 
S  J  C  J  and  S  • 

This  is  the  point  at  xifhich  the  investigations  of  this  report  stand. 
It  is  proved  directly,  from  an  inspection  of  the  infinite  determinants  G^,  S  , 

that 


(lU) 


f    .   2        \A               /_       .          \-     sinn  co 
(sm   n  co)  Z^        =     (2co  suinco  )   C     -s 


2  2 

(sin   nco)/^  -1  =  (cos    nco)  C     S 


(see  theorem  2  in  Section  II).  As  a  consequence,  four  analogues  of 
(7)  can  be  proved;  one  of  these,  for  instance,  is 

(15)  -  (2cosin  n  co)  C^  =  y^(^)  , 

(theorem  3,  Section  IV), 

Another  result  which  can  be  obtained  from  the  particular  nature  of 
the  determinants  C  ,  ,  .  .  ,  is  an  approach  to  the  inhomogeneous  Mathieu  equa- 
tion whose  right-hand  side  has  a  period  n  .  It  is  shown  in  Section  I  how  the 
method  of  infinitely  many  linear  equations  can  be  used  for  the  construction  of 
a  periodic  solution  of  the  inhomogeneous  equation  in  the  case  where  the  homo- 
geneous equation  does  not  have  such  a  solution. 

The  explicit  form  of  C  and  the  general  observations  of  Section  I 
about  the  determinants  of  this  type  reveal  that  the  left-hand  side  of  (13) 
is  an  analytic  function  of  to  which  can  be  written  as  a  sum  ha^/ing  the  form 

+00 

-2u)  Sinn  oo  +  >    s  

*—    n  u  +  n 


-  u  - 

Consequently  it  is  to  be  expected  that  -2w  sinnco  (C  -  1)  can  be  written  in 
the  form  of  a  finite  Fourier  integral,  since 

V2 

Because  of  (15)  one  may  expect  that  in  general,  for  any  value  of  x  and  any 
solution  of  y(x)  of  (l)  or  its  first  derivative,  a  similar  representation 
will  be  possible  if  certain  initial  terms  are  first  subtracted  from  y(x) . 
That  this  is  indeed  the  case  Is  stated  in  Section  V,  theorem  U.  This 
theorem  may  be  described  as  a  Fourier  representation  of  the  solutions  of 
(l)  with  respect  to  the  parameter  co. 

Finally,  in  Section  II  we  prove  the  well-known  theorem  about  the 
mutually  exclusive  nature  of  the  four  conditions  (10),  using  the  infinite 
determinants  C  ,  S  ,  ..«  .  This  theorem  is  the  negative  aspect  (in  the 
case  of  Mathieu's  equation)  of  the  problem  of  finding  equations  of  Hill's 
type  (8)  nhich  have  two  linearly  independent  solutions  of  the  same  period. 
In  section  VII,  the  connection  between  this  problem  and  the  problem  of  con- 
structing "transparent"  media  is  summarized  and  reviewed. 

Section  VI  indicates  the  possibilities  for  generalizing  the  results. 
The  results  obtained  for  Mathieu's  equation  could  be  proved  completely  for 
Hill's  equation  (8)  as  well  if  the  analogue  of  (lit)  could  be  proved.  Without 
that  only  some  results,  including  theorem  i;,  can  be  generalized. 

The  methods  used  in  this  report  are  based  upon  theorems  about  func- 
tions of  a  given  order  of  growth.  This  approach  to  the  theory  of  Hill's  equa- 
tion is  due  to  Schaefke. 

I  wish  to  thank  Professor  J.  Keller  and  Dr.  C.  J.  Bouwkamp  for  their 
valuable  help  in  preparing  this  report,  I  am  indepted  to  Professor  Keller  for 
the  method  of  Section  I  for  the  inversion  of  certain  infinite  matrices.  Dr. 
Bouwkamp  called  my  attention  to  several  papers  on  Mathieu  functions  and  Hill's 
equations  and  gave  me  his  unpublished  results  which  are  mentioned  in  Section 
VII, 
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I .  On  a  Special  Class  of  Infinite  Matrices  and  Determinants 

We  shall  consider  matrices  of  the  following  type:  The  elements  in 
the  main  diagonal  are  equal  to  unity,  and  the  only  other  elements  which  may 
be  different  from  zero  are  those  in  the  diagonals  adjacent  to  the  main 
diagonal.  The  most  general  matrix  of  this  type  can  be  indicated  by 


(1) 


1 

p. 

3 

0 

•         1 

0 

0 

• 

0 

0 

•        ^-3 

1 

P-2 

0 

0 

0 

0 

0 

q. 

2 

1 

p-1 

0 

0 

0 

0 

0 

'^-l 

1 

Po 

0 

0 

0 

0 

0 

% 

1 

Pi 

0 

0 

0 

0 

0 

% 

1 

.P2        • 

. 

. 

. 

• 

« 

..,  ■  *. , 

'•■              • 

We  will  denote  this  matrix  by  "  M(p,  .q^,)  and  its  determinant  (if  it  exists 

+  00     if       V 

in  the  sense  explained  below)  by  lrr)^(v^>%)}     the  subscript  +oo  and  the 
superscript  -oo  denote  the  range  of  rows  and  columns  or  the  range  of  values 
of  V  which  defines  the  set  of  parameters  p  ,q  , 

Similarly,  we  shall  denote  the  infinite  matrix 


(2) 


Po 

0 

0 

1 

Pi 

0 

^1 

1 

p 

by     ^K(p^,,a  ,)j       this  means  that  the  matrix  is  infinite  only  in  the  direction 

oo    V   V 

of  increasing  subscripts  denoting  roi^'s  and  columns. 
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If  the  rows  containing  p  ,p  ,p   ...  and  the  columns  containing 
^1  ^2   3 

Q,  jQ,  >Q,  >  •••  3^®  taken  out  of  M,  the  remaining  matrix  will  be  denoted  by 

\,\\     ^2  ^2  ^3  ^3 
and  the  corresponding  notation  will  be  used  for  determinants  and  their  sub- 
determinants. 

Definition  of  a  determinant  of  an  infinite  matrix.  Let  M  be  an  infin- 
ite matrix  and  let  A-i  ?  ^o>   ^■^'  *"  ^®  ^  sequence  of  finite  determinants 
formed  from  the  elements  of  M  (in  the  arrangement  in  which  they  occur  in  M) . 
These  detenninants  are  such  that 

(I)  The  main  diagonal  of  each  Aa  iZ   =  1,2,3,...)  is  part  of  the 
main  diagonal  of  M 
(II)  Every  element  of  M  appears  in  a  A.  and  every /\.  is  a  subdeter- 
minant  of  A/ .-i ,  and 

(III)  lim  A«  =  D  exists. 
/  ->oo  *» 

Then,  if  D  is  independent  of  the  particular  choice  of  the  sequence  >-A/,  we 

shall  call  D  the  value  of  the  determinant  of  M.  We  have: 

Lemma  1.  Let 

(3)  \   =  Pv  %  • 

Then  the  determinant  D  of  a  matrix  M  of  type  (l)  or  (2)  exists  if 

converges,  where  the  sum  is  taken  over  all  values  of  v  from  -oo  to  +oo  or 
(in  the  case  of  a   M)  from  0  to  oo.  The  value  of  D  depends  only  on  the  w^ 
and  not  on  the  P^,q^  themselves. 
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Ve  may  prove  this  lemma  by  observing  that  a  finite  subdeterminant. 


say  D,  can  be  shown  to  be 
•^  n 


O  "  *  XT * 

(ii)  D«l+V"w+yww       +x  www      +.., 

^      V        Z-         V^    Vg      ^^^^^       V,^    V^    V3 

where  the  si^ms  are  taken  over  v,v,,v„,v-,,    ...   and  wtiere  the  asterisk  jjnplies 

that  all  those  terms  which  involve  a  product  w  w     ^   of  two  consecutive  factors 
w     should  be  left  out  of  the  sum.     It  follows  that   the  infinite  product 

<5)  TT(i^  |wj) 

V 

Biajorizes  the  sum  of  the  absolute  values  of  the  terms  in  (U).     Then  the  con- 
vergence of  an  infinite  product  of  type  (5)  can  be  proven,  using  the  well- 
known  convergence  test,  and  this  proves  lemma  1, 

Using  the  formula  for  the  expansion  of  a   determinant  in  terms  of 
subdeterminants  belonging  to  the  elements  of  a  fixed  row,  we  derive  easily  the 
Reciirrence  relations; 

(7)  °D.     "t,         -^0     ''\.1'2,2--      '"     -»o'     ^. 

00  00        ^  00  '  '  00  00 

o 

(8)  D  »  lim       z 
^00  n 

n->oo 

where  z     is  defined  by  the  recurrence  relation 
n 

^^^     S+1  "  ^n  '  ^n  Vl   '   ^^'2  "  °'^-l  '  ^>   Zq  =  1>  ri-  0,1,2,...)    . 

Here  z^  _  is  actually  the  determinant  whose  elements  are  those  elements  in 
r-1 

_.D  which  appear  both  in  the  first  n  rows  and  in  the  first  n  columns, 
00 


-1  ° 

Now  we  shall  try  to  find  a  formal  inverse  M   of  a  matrix  M  =   M 

00 

of  type  (  2)  .  For  this  purpose  we  shall  use  a  method  found  by  Professor 
J.  Keller.  Let  Y,X  denote  the  infinite  matrices  indicated  by 

1   X    0   0... 


(10) 


1 

0 

0 

0 

^0 

1 

0 

0 

0 

^1 

1 

0 

0 

0 

^2 

1 

Y  , 


0 

1 

^ 

0 

0 

0 

1 

X 

0 

0 

0 

1 

2  • 


X. 


Then  a  formal  multiplication  of  these  matrices  gives 


0    0 


(11)  TL 


y  1+y  X   X, 
''o   •'o  0   i. 

0     7^  ■'■^yi^ 


0    0 


0    0 


0    0 


72  1*72^  ^3 


73  ^^yf-2 


'o>'=l>^2- 
in  the  main  diagonal,  and  postulating  that  EYX  shall  be  equal  to  M,  we  find 

the  relations 


e     "I 
o 


(12) 


'1^0  =   %   »     ^l^^^Vo^  =  ^>     ®1^  "  Pi 


%.l  V   %  '  ^.1^^  *  ^n^n^   =  ^'   ^n+1  Vl  "  Pn.l 
We  derive  from  (12)   : 

^^^^  ^n  -  %/«n.l     '  Vl  "  Pn.l/^n.l  * 


(lU) 


'n+1 


1  -  w  /e 


1  -  X 


h^ 
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Since  e     »  1  and  e,   =  1-w    we  find  for  e  ^^  the  continued  fraction 


(i5>  vi"i-  r^ 


1  -  w  , 

n-1 


1  -  V2 


^1 


"-'% 


From  (15)  we  see  that  y  and  x  ^-  are  determined  completely  by  (13),  provided 

that  none  of  the  e  equals  aero.  We  can  now  find  a  formal  inverse  of  M  by 
n 

observing  that 

^  "^0  '  ^o\  '   "^o^^2  '     \   ^  ^2^3  >    •    '    '  .    •    •    ■ 

0    1      —X,   J   '*^  ^T  ^O    J  "  "^  ^'P  ^"^      ^  .  .   . 

(16)  X"-^  =     °   °     ■••     "^2        ^2^3     ^  •  •  * 
0   0     0      1         -x^        ... 

0   0     0      0  1       ... 

We  simmarize  the  result  by  stating 

Lemma  2.  Let  e  be  defined  by  (l5)  or  by  e^  =  ^^/"^^^i    whsre  z^^  is 

defined  by  (9).  Assume  that  lim  z  exists  and  is  different  from  zero.  Let 

n->oo 

W  WW,  WW-W^o 

/,r,\      „         -,  n  n     n+l  n     n+1     n+2     ^ 

(1?)   ^n '  1  ^  rr—  *  — 2 "  — r-~5 -  -  -  > 

"  "*^    ^n^n+l^n+2    «n%*l^n+2  ^n+3 
or,  recursively,  let 

(C  -1)  Z  7   n 

^  n  '  n^n+l 


(18) 


n+1    w  2  T  z 
n  ii-1  n 


For  n,m  =  0,1,2,  .  .  .  let 

(19)  B„=  (q^q^  ...qJ/Vl>  ^.i  «  0  • 
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Finally,  let  us  assume  that  all  of  the  e  are  different  from  zero.  Now, 
consider  the  matrix  M  given  by  (2),  i.e., 

(20)  M  =  i[i^^)  (n,ra  -  0,1,2,  ...) 
The  inverse  of  this  matrix  can  be  written  in  the  form 

(21)  M"^  =  (ti*   )  , 
where 

(22)  M*   -  C  /e 

^  n,ra    n'  n 


and,  forX  =  1,2,  ...  , 

(23) 

*  i 

HA-  (-1)   C  6      B   fl  T  /  (B  T  e  )  . 
■^n+Jt  ,n   ^  '    n+*   n+i  -1  '   n-1  n' 

If  one  of  the  e  (or  z  )  vanishes,  these  formulas  become  meaning- 
n     n         ' 

less  and  have  to  be  replaced  by  more  complicated  ones.  Since  it  is  easy  to 

show  that  the  sum  of  the  squares  of  the  modioli  of  the  elements  in  X~  -  I  and 

Y~  -  I  converge  (I  denotes  the  identity),  we  can  conclude  that  M   is  a 

bounded  matrix  if  lira  z  exists  and  is  different  from  zero  and  if  5~  |w  I  is 
n  *^ —   n  ' 

n->oo 

finite . 

We  also  mention  the  following  fact,  >iiich  is  proved  easily: 

Lemma  3 »  If  in  the  matrix  M  one  element  in  the  main  diagonal,  say 
(jL   ,  is  replaced  by  zero,  M"  still  exists  and  can  be  expressed  in  the 

manner  of  Lemma  2  if  we  determine  e,  ,x,  ,y,  ,  by  the  relations 

^^^^  \   =  ^-i/^k-i " "  vy^k-i 

\     -  Pk/^k"  -Pk«k-i/Vi, 

Here  we  assume  that  w,  -,  •=  p,  ,  q,  ,  j^  0 ,   and  that  the  assumptions  stated  in 
lemma  2  still  hold.  For  all  values  of  nji  k,   the  y^^.^^,  e^,  x^  have  to  be 
determined  as  before. 
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We  can  use  the  results  of  this  section  to  find  a  periodic  solu- 
tion of  the  inhoinogeneous  Mathieu  equation 

(25)       y"  +  h{(/   +  2t  cos  2x)y=  f(x)  , 

where  f(x)  is,  for  instance,  an  even  fxmction  having  period  n  and  having  a 
continuous  first  derivative.  We  shall  expect  (25)  to  have  a  uniquely  deter- 


mined even  solution  z(x)  whose  period  is  n,  if  the  homogeneous  equation 

(25)  y"  -f  h{(i>^   +  2t  cos  2x)  y  =  0 

does  not  have  such  a  solution.  We  can  find  2(x)  by  putting 


00 

(27)  z(x)  "    y      £  cos  2nx 

n=0   " 

00 

(28)  f(x)  "    ^      b  cos  2nx  . 

n=0    " 

Then  we  find  the  follovang  equations  for  the  4  * 

-  co^4   -  t  Cn  «  -  b 

^o      ^1  0 

(29)  -  2t  ?^  +  {1-J)K-^-  t  $2   =  '\  1^ 

2 
-  tC;L  /2^  +  (1  -  ^  )  ^2  -  -t^  /2^   -  -  ^2  /^^ 

The  determinant  of  t  he  linear  equations  is  exactly 

(sinnw)  C     , 

IT  + 

where  C^  is  defined  by  (II, 20).  The  vanishing  of  oo  sinn  w  C  is  a  necessary 

and  sufficient  condition  for  the  existence  of  an  even  solution  of  period  n 
of  equation  (2.b).  If  no  such  solution  exists,  we  can  apply  the  results  of 
this  section  to  the  system  (29),  and  if  none  of  the  quantities  e  vanishes, 
it  is  not  difficult  to  verify  that  the  solution  of  (29)  for  the  ^  is  such 
that  z(x)  has  a  continuous  second  derivative. 
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II.  Infinite  Determinants  Related  to  Periodic 
Solutions  of  Mathieu's  Equation 

We  shall  write  Mathieu's  equation  in  the  form 

(1)  y"  +  k{oi^   +  2t  cos  2x)  y  «  0 

and  we  shall  state  the  conditions  under  which,  for  a  given  a,  yexp(-2iax)  is 
a  periodic  function  having  period  n   for  at  least  one  solution  y  of  ( 1  ) .  The 
customary  form  of  (l)  is 

(2)  y  +  (a  -  2qcos  2x)  y  «  0, 
and  2a  is  usually  denoted  by  ii. 

If  o  =  m/n,  wJ^re  m,n  aire  co-prime  integers,  and  if  yexp(-2ittx)  is 
periodic  and  has  the  period  n,  then  y  itself  is  periodic  and  has  the  period 
nn  e  It  is  known  (and  will  be  proved  again  here)  that  for  n  -  1,2,  there  do 
not  exist  two  linearly  independent  solutions  of  (l)  both  of  which  have  the 
period  nv..     If  n  >  2,  however,  there  will  be  two  such  solutions  for  appro- 
priately chosen  values  of  w  and  t. 

It  is  well-known  that  for  given  values  of  co,  t  the  value  of  a  can 
be  computed  from  either  (U)  or  (15). 

Let  Yt ^yp  denote  the  solutions  of  (1)  which  satisfy  the  initial 

conditions 


(3) 


Then 


y-^Co)  -  1,  y^^Co)  -  0 

ygCo)  »  0,         ygCo)  -  1. 


(U)  cos2na  •  y,(jijco,  t)  , 

where  we  mention  the  parameters  oo,  t  explicitly  in  the  function  y,  (xjco,t)  . 

Since  exp2Tiia,  rather  than  a   itself,  describes  the  behavior  of  y,  we  may 

restrict  a  in  such  a  way  that 

(5)  0  ^  Rea<  1, 
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Then  (U)  gives  two  possible  values  for  a  if  the  right-hand  side  of  the  equa- 
tion is  known,  since  if  a  is  a  solution  of  (U)  satisfying  (5),  then  1  -  a  is 
also  such  a  solution.     However,  if  o  «=  ^  or  c  «=  0  there  will  be  only  one 

value  of  a  satisfying  (U)  and  (5).     It  is  known  that  for  every  a  satisfying 
(h) ,{S) }  there  actually  exists  a  solution  y  of  (l)   such  that 
y(x+r)   =  exp(2nai)  y(x).     We  shall  prove  now: 

LemiTia  U.  Let  a  be  a  solution  of  (U),   (5)   and  let  a  be  different  from  zero 
and  from  1/2.     Let  y,2  be  solutions  of  (l)    such  that 

/         N  2Tiia     /    V        /         N         2ni(l-G)      /    ^ 

y(x+Ti)     »  e  y(x),   2(x+ii)  =  e       ^  z(x)    . 

Then  y(x)   and  z(x)   are  lineai^ly  dependent.     Proof:     If  y  and  z  are  linearly 
dependent  we  ntay  assume  that  y  "  p  z  ,  where     p  j<  0  is  a  constant .     This  gives 

/         s         2nia     /    X  ,        ->  -2nia     ,■    v 

yCx+n)  =  e  y(x)  =     PzCx+rJ  »     p  e  a(x;. 


Therefore 


y(x)   =     pz(x)  »  p  e"  z(x) 


This  iiriplies  that  exp  (linia)   »=  1  and  therefore  that  either  a  «  0  or  a  »  -^  , 
which  contradicts  our   assumption. 

As  an  immediate  consequence  of  lemma  k  we  have: 
Theorem  1 .     If  Mathieu's  equation   (l)  has  a  periodic  solution  with  a 
period  nn,  n  >  2,   all  of  its  solutions  have  this  period. 

The  possibility  that  all  solutions  of  (1)  may  be  periodic  has 

been  mentioned  by  Hille  [2^    .     for  even  values  of  n,   the  existence  of  two 

solutions  of  period  n  for  appropriate  values  of  co  and  t  was  proved  by  Foole 

L?]     .     Since  it  is  known'    that  (U)  can  be  satisfied  for  any  real  value  of  a 

and  infinitely  many  pairs  cf  values  of  co  and  t.   Theorem  I  provides  examples 

for   "transparent  media"  ,  i,e.,  for  the  Ramsauer  effect,   in  which  a  plane 

wave  passes  through  a  finite  layer  of  a  refracting  medium  without  reflection. 

*  y,TT'>w,t)  =  cos  2na  is  a  continuous  function  of  co  and  t.     This  follows 

directly  from  Picard's  method  of  solving  a  differential  equations  by  itera- 
tion.    The  general  theory  of  Mathieu's  equation  shows  (see  e.g.  McLachlan 
[5]that  y,  (Tt,u3,t)  is  real  for  real  values  of  t  and  co  and  that  for  a  fixed 

t  it  takes  on  the  values  +1  and  -1  infinitely  many  times  for  real  values 
of  ffi.  Therefore,  cos2na  taJkee  on  energy  values  between  +1  and  -1  infin- 
itely many  times. 
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For  a  more  detailed  account  of  tMs  matter  see  Section  VII, 

The  following  statements  about  solutions  of  (l)  which  are  periodic 
and  have  the  period  n  or  2n  are  standard  results  in  the  theory  of  Mathieu  func- 
tions. 

I.  A  solution  which  has  the  period  n  is  either  even  or  odd  and  can  be 
represented  in  one  of  the  following  two  ways: 

00  00 

(6)  ce  (x)  =  y"  A^„  cos  2rx  ,    se(x)   =  T"  B^  sin  2rx  , 

r=0  '^^  T=l    ^^ 

II.  A  solution  which  has  the  period  2n  but  not  the  period  n  is  either 
even  or  odd  and  can  be  represented  in  one  of  the  following  two  ways: 

00  00 

(7)  ce__(x)  «  J3  k^^^^   cos  (2r+l)x  ,  se_(x)  ^  ^^  B^    sin  (2r+l)  x. 

r=0  ~     r=0 

III.  If  for  a  certain  pair  of  values  of  ai  and  t  there  exists  a  periodic 
solution  of  (1)  which  is  of  type  (6)  or  (7),  the  multiples  of  this  solution  are 
the  only  periodic  solutions  of  (l). 

We  shall  now  prove  I,  II  and  III  by  using  infinite  determinants. 
Assuming  that  (1)  has  a  solution  of  the  type 

+00 


(8)       y  =  e^^'^  YL     '     '^^'' 


c  e 

m 

m=-oo 
to  which  two  term-by-term  differentiations  can  be  applied,   the  c     have  to 
satisfy  the  equations 

2   /         x2 


(9^     ^  ^  Vl  *  ^C"  -(m+a)    ]Cj^  +  U  +  litc  «     0. 


m+1 
(m  =  0,  :r  1,  !  2,   ...   ) 

(.• 
2      2 
Dividing  (9)  by  U(a)  -m  )  ,   the  determinant  of  the  resulting  system  of  linear 

equations 

2  2 

(^n\  t  oi  -    (m-a)  t 

^^°^     -T-2     Vl   "     — —     °m   '     -2 2     °m.l     '  ^ 

CO   -m  CO     -  m  co     -  m 

exists  and  shall  be  denoted  by 
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We  write 

(12)  Aq  (",t)   =  A(05  ^,t)  , 

and  then  there  exists  the  relation  (see  for  instance  Whittaker  and  Watson jlU"}  ) 

(13)  A(a=,  co,t)   -   A>,t)   -    (fl^)'      . 

Here  A  (co,t)  can  be  written  as  a  determinant  "oo^^Pv»'^v^  °"^  ^^®  (l«l)» 

(li+)        Pv  "  %..!   =  -i~2     '      ■ 

m  -0) 

A  necessary  and  sufficient  condition  for  the  existence  of  a  solution  of  (l) 
of  type  (8)  which  can  be  differentiated  twice  terra  by  terra  is  the  vanishing 

of  sin  no)  A(a,>jo,'t  )  .  Therefore  we  obtain  from  (u)  and  (13): 


(15)   cos2Tia»  y-,(n,w,t)  »  1  -  2  A^(a),t)  sin 


2 


If  (l)  has  a  periodic  solution  ce  (x)  of  the  type  (6),  the  A-  must 
satisfy  the  set  of  linear  equations       ^  _ 

Uco^  A^  +  Ut  Aj  «  0 
...v         Qt  A^  *   U(co^-l)  A^  +  Ut  A^  =  0 

Ut  A^^_^   +  U(oo^-r^)  Aj^  +  Ut  A2j.^2  "  °  * 

(r  -  2,3,U,  ...)  . 

Similarly,  the  existence  of  periodic  solutions  of  the  type  se^,ce_,   se_  leads 
to  the  following  systems  of  linear  equations: 

U(co^-l)  Bg  +  UtB^  «=  0, 

(17)  P     p 

Ut  B2^_2  *  i*(^  ~r   )   ^2^  *  Ut,  B2^^2  '  ^    > 

(r  =  2,3, U,    ...)i 
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(18)  U  {["^  -  (|)^  J  +  tj  A^  +  Ut  A3  -  0  , 

2  1  2 

^t  A^^.i  +  U    [co     -(r  +  ^)     J   A^^^^  +  at  A2^^3 

(r  =  1,2,3,...); 


0   , 


(19) 


k  Ifco'-Ci)'  ]   »  t}  B^  +  Ut  B3  -  0   , 
2  1  ^ 

(r  -  1,2,3,...)    . 


0    , 


If  we  cLLvide  the  equations  (16),  (17),  (I8),  (19)  by  certain  factors  in  such 
a  manner  that  the  matrix  of  the  resulting  system  of  linear  equations  becomes 
the  unit  matrix  for  t  =  0,  we  arrive  at  the  following  infinite  determinsurits 


(20) 


C.    = 


1 
2t 

J2 


t/co 

1 


2    „2 
CO   -2 


0 

t 
co^l 


0    , 
0    , 


2  „2 
CO  -2 


2     ^2 
CO     -3 


(21) 


> 

CO  -1 

^ 

> 

^   » 

t 

> 

1 

» 

t 

> 

0    , 

• 

2    ^2 

CO   -2 

0 

t 

CO   -3 

» 

1 

y 

t 

CO  -3 

> 

•   • 

0 

0 

t 

) 

1 

> 

co^-U" 

(22)  C  = 


t 


co2.(|)' 


-  17 

__t 


o2-(|)^ 


0,0, 


^^-(l)' 


,  0 


.^-(l) 


.^-(l) 


(23)  b 


1  - 


2  rl^' 


»^-(|)' 


2  An 

J  -(^) 


co2-(|)' 


,    0   , 


.^-(|) 


T  ' 


.2.(1)' 


o^^-il) 


Of  these  determinants,  C  and  S  are  of  the  type    D  of  Section  I,  and  C  and 
S  differ  from  this  type  only  in  the  first  element  of  the  main  diagonal. 
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We  shall  now  prove: 
Theorem  2. 

(2U)      sin^no)  A^C^jt)  •  (sinncoC^  )  (sinncoS^  ) 

(25)  sin   nco  £^  (co,t)  J   -  1  =  sin    nco  {"AC^jw,  t)3   =  -  cos    nco  C_  S_    . 

If  t  is  a  fixed  ntomber  different  from  zero,  no  more  than  one 
of  the  functions 

(sinno))  C  ,  (sinnw)  S  ,  (cos nco)  C  (cosnci))  S 

vanish  for  any  given  value  of  co.  Equation  (2U)  establishes  the  well-known 
fact  that  the  condition  sin  noi  ^  (a),t)  »  0  for  the  existence  of  a  solution 

of  period  n  is  satisfied  if  and  only  if  one  of  the  conditions  (sinnw)C  =0  or 
(sinn(a)S  «  0  is  satisfied,  that  is,  if  a  solution  ce  or  se^  exists.  Equa- 
tion (25)  establishes  the  corresponding  fact  for  the  solutions  of  period  2n. 
However,  (2U)  and  (25)  state  more  than  this:  they  establish  a  factorization  of 
an  entire  function  (namely  sin  nco  Z^  (oi),t)  )  into  the  product  of  two  other 
entire  functions.  We  can  derive  (22)  from  the  recurrence  formulas  (1.6),  (1.7) « 
We  shall  use  the  notations  of  Section  I  and  in  particular  we  shall  denote  sub- 
determinants  by  double  subscripts.  Then,  if  we  put 

(26)  P_^  -  P^  =  t/(co^.v2)  ,  q^  »  p^^^  ,  w^  =  p^q^  , 

remembering  that  p=p,w=w,,we  find  that  the  subdeterminant  of 

^  (co,t)  which  corresponds  to   ~,  D  in  (1.6)  would  be  C  if  in  the  first 

2  2 

column  of  C^  we  substitute  t/(co  -1)  for  2t(co  -1).  Since  (1.7)  gives 

C  =  S  -  2w  S  ,  T 

(27)  +    +     0  +,1,1 

S    "   S    T   T   -  W,   S     ,    ,   o   O 

+    +,1,1    1  +,1,1;2,2 
we  find  from  (1.6)  that 

(28)  A^io^.t)   =  (S^.  w^S^^^^)  (S^^^^-w^  S^^^^^^^)   -  w^S^S^^^^ 
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which  proves  (2U).  To  prove  (25)  we  observe  first  that 

1~ 


f  .    »  2      +00    2    2 

(29)     iii.^1   .  -  rr  '^  -^ 


n»-co   2  /„^\ 
CO  ~  (n+^; 


if  the  product  is  defined  as 

lim 
N->co 


+N 


n  . 

n=-N 


-1       2   2    2     1 
Multiplying  equation  (10)  (for  a  «  j)  by  (co  -m  )/(w  -(m+j)  ),  we  find  from 

(29)  that 

(30)  ^   -  tq^nca  /i!^|j  «,t) 

is  an  infinite  determinant  of  type   (l.l)   if  we  specify 

2      •      -,     2 

(31)  P^  -  t/|^co  -  (v+^  )     J   »     %  "  Pv-1  '  ''v  "  Pv%  • 

We  then  have 

(32)  p_^  -  P^.i  ,   q_y  -  Py.2»     w_i  -  Po      • 
Defining  the  determinant  T  by 

(33)  T  - 


1        P, 
Pi 


o 
1  P 


0  0  0 

0  0 


0 


P2       0 


where  the  p  are  given  by  (31),  we  find  that 


(3U)      S_  .  T  -  p^  T^^^  , 


C  -  T  +  p 


o  ^1,1  • 
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-c» 
On  the  other  hand,  if  we  substitute  /^   for     D  in  (1.6')  we  obtain 


-00  -00  1  2 

(35)  .2  D  «  T  ,     -3  ^  •=  ^1,1  '  cx.''^  =  ^1,1  '   -  ^  -  ^l,l;2,2 

and  therefore,   from  (1.6)^   (!.?)>  and  (32)    : 
*  2     2 

(36)  ^  =^  V^o  V-^^  V;^,2 


^  -  Po  '^1,1     =  S.  C_   . 


This  proves  equation  (?5)  of  Theorem  2. 

To  prove  the  statement  of  Theorem  2  about  the  non-existence  of 
common  zeros  of  the  four  functions  (sin  nco)C     ,  etc.,  we  observe  first 
that  (2U)  and  (25)  exclude  the  possibility  of  (sin  Tia>)C     or  ( sinner) C 

having  a  zero  in  coiranon  vdth  (cos  nco)C     or  (cos  nai)S     ,     The  fact  that 

S_^  and  C^  cannot  have  a  zero  in  comnon  unless  t  =  0  can  be  derived  from  (2?) 

as  follows:     If  S^  «  C^  =  0  and  t  ^  0,  we  would  also  have  w     f*  0,  w,   ^  0  and 

consequently  S     -.^-S     -,   t   r,  r.  '  0 .     Byan  application  of  (1.7)   to  S     .,   -, 

and  S  -,1,152,2  this  would  lead  to  S^.-,   ^  .„„.-,,»  0  and  so  on,  which  is  im- 
+  ■•■ji-,Xj<;,<i,^,j> 

possible  since  the  subdeterminants  of  this  type  approach  unity,  llie  possi- 
bility that  (sinnco)  C  and  (sinitoi)  S^  could  have  a  common  root  at  w  =  mn  , 
m  ■  0,i  l,i2,.  .  can  be  excluded  for,  if  there  were  such  a  common  root,  the 
row  of  (sinnco)  C     which  is  labeled  by   |mj  would  become,  for  co-o  i  m  (m/O)^ 

0,  0,     ij^IL  ,     0,     itr^    ,^  ,    _    _ 
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where  the  zero  in  the  center  is  the  element  of  the  main  diagonal.  The 
corresponding  recurrence  formula  would  then  read 

^+,1^;  ...  i  |m|-2,  |m|-2  "  ^^+,1,1^..-;  \M,   1^1, 

where  X  ^  0.     The  result  is  a  contradiction  as  before. 

The  non-existence  of  common  aeros  of  (costio))  C_  and  (cosnco)  S^  can 
be  proved  in  the  same  manner. 

III.  A  Majorant  for  Mathieu's  Determinants 

There  are  two  methods  of  obtaining  a  majorant  for  the  absolute 
value  of  a  determinant  of  type  (l,l)  or  (1.2).     The  first  is  to  apply 
(1.5).  which  givesTI(l+    |w   |)   as  an  upper  bound  for  the  detei'irdnant.     The 
second  method  is  to  apply  Hadamard's  inequality  which  states  that  the 
modulus  of  ttie  determinant  of  the  elements     a  is  at  most  equal  to  the 

product  of  all  s   .  where 


Also,  we  may  observe  that  for  a  determinant  D  of  type   (l.l)  or  (1.2), 

>  Id-  i| 


(2)  TT  (1  *  \\\)-  1 


since  the  terms  on  the  left-hand  side  in  (2)  majorize  the  moduli  of  the 
terms  of  D-1  if  D  is  expanded  in  a  sum  of  type  (I.U).  Applying  (1.5)  to 
A  (w,t),  which  is  defined  by  (Jli.l2) ,(J.lh) ,   we  find 


.2        ,  <   1.    2 


00 


(3)    \^  (co,  t)sin  nw   j   =    [sin  nco    |  "rf        1  + 

n=-oo 


2 


t 
2    2x,,   .,x2  .2 


(n  -co  )  |(n+l)  -03      I 


r ^ 


■fr       /     i(n^-co')(n^l)^-co^)f*|tr 
n=-oo       I  n         (n+l) 
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where  n     ■  n  for  n  ?<  0,  0     «  1,  and  (n+1)     »  n+1  for  n  j<  -  1,  0     -  i.     Now 

(U)  I   (n^-co^)   ((n+D^-co^)   |   +    Itj^  <  (n^4   jc^^  |  ^    |tW  |)((n+l)^  +  |oF  | )  , 

and  therefore  we  find  from  (3)    : 


<     Isinhw    U  |<o   I  •>•    It  /ml  J       slnhtilcoj  J 


(5)  I A     (co,t)   sin'^no)  |   .     _— ^ 

I  sin  na>| 

Combining  (5)  and  (2)  we  find  that,   for  a  fixed  value  of  t  and  for  co  «  iX, 
X  real, 

(6)  lljn       |A   (iX,t)    1-1. 
X->  00    '  ' 

To  apply  Hadamard's  inequality  we  observe  first  that  the  expansion 

of  sin  n CO  into  an  infinite  product  shows  that  the  n   row  (n  ■  0,11,  12  ,  ...) 

A         2 
ofZa  (a),t)sin  tico  can  be  written  as  follows: 

For  n  y'O  , 

?  2 

.    .    .   ,  0,  0,  0^-t/n     ,  1  -  ^  ,  -  t/n^   ,0,0,0,    .    .   ,,and  for  n  -  0^ 

n 

2  2  2         2 

.    .    .    ,  0^  0  ^  On  t   ,   n  CO     ,  n  t,  0,  0,  0,    ,    .    . 

Therefore  we  find  from  (2): 


(7)    |A  (co,t)sin^nco  |  <  n^(  |co|^  +  21t|) 


{aC'-'^)'} 


1 

sinh^  Ti   r  |w|^  +  2|t|  3^     . 
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Applying  the  same  method  to  C  ,  S  ,  C  ,  S_  of  Section  II,  we  obtain  results 
which  may  be  6\immarized  as  follows: 

Lemma  5«  For  all  values  of  co  and  t  we  have  the  inequalities 

1 


(8)      |sinnco|   |C 


J   <  M^2Lti!.,^,[;,(H2,,,,,)l  ]  ^ 


(9)      jSinjw_  g 


<     arJiCn  (|co^|   *  2|t|)^J 

1  ' 

i\^^\   +  2|t|  f 

1 

(10)  icosnco  ||C_|     t  cosh  [ti(  |ai|^  +  2!t|)^J    , 

1 

(11)  |cosna)|jS_l       <  cosh  [ji(  ico|^  +2  jt  |)^  J    . 

The  inequality  (7)  is  an  immediate  consequence  of  (8)  and  (9).     Lemma  5, 
together  with  the  relation  (6),  implies  a  result  foiind  by  Sohaefke 
(theorem  VII,  p.  39  in  |llj  )  about  the  order  of  growth  of  (sin  nw)  /^  («t)  . 

IV.     Expressions  in  Terms  of  Infinite  Determinants  for 
the  Values  of  the  Solutions  of  Mathieu's  Eqxiation  for 

X  ■  »  and  X  "  n   . 


In  this  section  the  following  result  will  be  proved: 

Theorem  3  •  Let  y,(x5a5,t)   and  yp(x;m,t)  denote  the  even  and 

odd  solutions  of  Mathieu's  equations  (ll.l)   and  assume  with  the  initial 
values  (1,0)   and  (0,1)  for  y,y'   at  x  »  0.     Let  G^,  S^,  C_,   S_  be  the  in- 
finite determinants  defined  by  (II. 20)  to   (11.23).     Then 

(1)  ^l^I'  ^»^^     "  ^°°®  ""^  ^- 

(2)  y]_(^>'^>t)      »  -2w(stnnco)  C^ 
f-,\  /  n       .  s  sin  no)    „ 

(3)  ygCji^^^t)     -  -^j;^-  s^ 

(U)  y2  (^j"  >t)     -  (cosnoj)  S_ 
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These   formulas  may  be  supplemented  by  the  elementary  relations 

(5)  2y^(|}  co,t)  y^Cp  <^,t)     -  y-,(Tii«,t)     -    1 

(6)  2y^(f}co,t)  yj  (f;a),t)       =  y^(n;<o,t)     +     1 

(7)  7jin<^»'^)  '  y2^^><^>^'>  » 

of  which  (5)   and  (6)  have  been  found  by  Schaefke  [[loj     and  (7)   is  stated 
in  McLachlan's  book   l5^ . 

To  prove  (l)  to  (U)  we  must  prove  the  foiled ng  facts; 

1.  Necessary  and  sufficient  conditions  for  Mathieu's  equation  (ll.l) 
to  have  solutions  of  the  following  types: 

even,  with  period  2n 
even,  iri. th  period  n 
odd,  with  period  n 
odd,  with  period  2n 


(8) 


respectively. 

Proof:     Assume  for  instance   that  y(.x)   «  y(x+  2n)   and  y(-x)   ■=  y(x) . 

Then  we  derive  from  y  «  c,y,(x)   <■  Cpy^Cx)  (where  c,  ,0^  are  constant 
values)   that  c,   =  y (O)  ,  Cp  ■  0   . 

Therefore  y  -  y.,(x)   is  the  even  solution  with  period  2n.     Because 
of  the  periodicity  of  the  coefficients  of  Mathieu's  equation,  y(x+n)  is  a3.so 
a  solution  if  y(x)  is  one. 

Therefore  we  find  that 

yj^(x+n)   -  y^(a)  y-j^(x)  ♦  y^(n)  Jg^^^ 

y2(x+T[)  -  ygCn)  J-j,^x.)  *  ygCTi)  y^Cx)    . 


yi(f;co,t) 

=     0 

j^{pui,t) 

»     0 

72(^5'^,*) 

-     0 

j^^r^^.t) 

«     0 
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The  matrix 


/yi(n)   yi(")'\ 
\y2(^)    y2(")/ 


has  a  determinant  which  is  the  Wronskian  of  the  fundamental  solutions  y,(x) 

V  2 

ajid  yp(x)  and  is  equal  to  unity,     M     must  be  a  matrix  in  which  the  elenenta 

of  the  first  row  are  1  and  0  since  y^(x+2n)  »  y^(x) ,     From  these  conditions 

we   can  derive  the  fact  that 

7-iM    \jAii)  +  yJn)    I   -  2,  y  (n)  »  0,     y  (n)         »     1  . 
112  1  I  1  , 

Therefore  we  have  y,(Tt)  =  -1,  y^C")  "  -!•  (We  discard  the  other  possibility, 
y-i(^)  =  yp('i)  ■  Ijbecause  it  would  mean  that  y.,  (x+it)  =  y,  (x)  and  hence  that 
y,  would  be  already  of  period  n.)  Therefore  we  find 

y^(x+n)  »  -y-^ix)   . 
For  X  "  -Tt/2,  we  find  that  y,(-n/2)  =  y,(Ti/2)  =  -y  (ii/2)  =  0.  On  the  other 
hand  if  y,  (ti/2)  =  0,  we  also  have  y-(-n/2)  =  0  ,  and  since  there  must  exist 
a  relation 

y^(x+n)  "  r^y-^M   +  r272^^^   > 

we  find,  putting  x  «  0  and  x  »  -  ■^  ,   that 

Yi  -  y^C")   »  Yg  =  0   • 

Therefore 

Jj^iO)  =  y-j^(-n+n)  »  y^(n)  y-j^(-r)  -|y-j^(Ti)^         «     1  , 

and  therefore  y,(x+ii)   =  -y(x)  unless  Yt   "  1»  iJi  which  case  y,(x)  has  the 
period  n.     But  then 

y^(x+n)  -  y^(x) 

gives  for  X  =  -  •^ 

y^(Tt/2)   -  y^(-  f)    . 
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Since  y, (x)  is  an  odd  function  of  x.  we  would  have  y, (n/2)  -  0^ which  is 
incompatible  with  y-,(n/2)  «  0.  Therefore,  y^(x+n)  -  -y-,(x)  follows  from 
y,  (n/2)  "0,  and  this  completes  the  proof  of  our  statement  1, 

The  next  statement  which  we  shall  prove  is 

2 
2.  Let  X  »  ItciJ  and  consider  the  equations  (8),  for  a  fixed  real 

value  of  t,  as  equations  for  X.  Then  all  the  roots  of  these  equations  are 

siiaple  ones. 

Proof;  We  use  the  fact  that  for  y-,(Ti/2)  =  0, 

^1  (^'/^^   o   ^ 

Equation  (9)   and  similar  formulas  for  the  other  functions  in  (8)  can  easily 
be  derived  from  Mathieu's  equation.     Now  if  t  is  real,  (  y-i(C)J       is  non- 
negative,   and  this  proves  statement  2. 

Finally,  we  prove  the  following  statement  about  the  left-hand 
sides  in  (8); 

3»     Let  X    be  any  fixed  (real  or  complex)  value  of  x  ard  let 
y(x;cD,t)  be  any  solution  of  Mathieu's  equation  (II, l).     Then  there  exist 
nvimbers  A  and  B  independent  of  co  such  that 

(10)  |y(x^;co,t)    t  <  A  e^l^l    . 

The  proof  is  part  of  the  wsll-known  method  of  solving  (ll.l)  by 
Iteration.     Choosing  for  ^  (x)   a  solution  of  y"  +  Ijca  y  =  0  which  satis- 
fies the  correct  initial  conditions  and  defining  y  (x)  for  n  »  1,2,3,    ... 
as  that  solution  of 

which  satisfies  the  correct  initial  conditions,  we  can  easily  see  that 

y  =  lim  iW        satisfies   (lO). 
n->oo  ^ 
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The  corresponding  facts  concerning  (cosnoj)  C  etc  are: 

la.  A  necessary  and  sufficient  condition  for  Mathieu's  equation 
to  have  an  even  solution  of  period  2n  (but  not  of  period  n)   is 

(lOa)  (cosjico)  C  »  0  . 

2a.  For  any  fixed  real  value  of  t,  the  roots  of  (IDai)  are  simple 

2 
ones,  if  o)  is  introduced  as  a  new  variable. 

3a.   (cosHoo  )  C_  satisfies  (III. 10), 

Statement  la.  has  been  proved  elsewhere;  see  for  instance 
Whittaker  and  Watson  ^lU J  .  Statement  3a.  was  proved  in  Section  III. 
Statement  2a.  can  be  proved  as  follows: 

Using  (11.15),  (11.25)  and  (6)  we  find: 
(11)      yi(Ttjco,t)  +  1  =  2  -  2  A^(co,t)  sin  nco 

'  2   cos  TicD  S_C_  -  2y^(^;co,t)  y2(^jco,t)  . 

Now  y^(r/2)  vanishes  if  and  only  if  (cosnco)  C^  vanishes,   and  yp(n/2)  vanishes 

if  and  only  if  (cos  nco)  S_  vanishes,   since  the  vanishing  of  any  of  these  quantities 
is  a  necessary  and  sufficient  condition  for  the  existence  of  a  solution  of 
Mathieu's  equation  with  one  of  the  properties  described  before  eq.   (8).     Now 
y,(n/2;a>,t)   and  y„(n/2;co,t)  have  simple  zeros  (with  respect  to  oj  )  for  any 

given  real  value  of  t.     Therefore  the  quotients 
,^-s  (cosnco)  C  (cosTtco)  S 


y^(n/2;co,t)  y2(V2>o>,t) 


2 

are  entire  functions  of  co  .     Since  their  product  is  unity,  neither  of  them 

can  have  any  zeros.     Therefore,  the  multiplicity  of  the  zeros  of  the  numerators 
cannot  be   higher  than  the  multiplicity  of  the  zeros  of  the    denominators.     This 
proves  statement  2a. 
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To  prove  equations  (1)  to  (U)  of  theorem  3}  we  csi  apply  a  theorem 
in  the  theory  of  entire  (or  meromorphic )  functions  the  proof  of  which  is  found 
for  instance  in  Nevanlinna  ["tJ  (p. 205-213)  or  Titchmarsh  ^13l  .  Kiis  theorem 
shows  that  from  statements  3  and  3a  the  following  representations  for 
y-(w/2jco,t)  and  (cos  nw)  C_  can  be  derived: 

2 
(13)  yi(V25co,t)  -  a)2«a(t)  IT  (l  -  f-  ) 

^  V  =1     °v 

(lU)  (costico)  C  -  (o^^'a'Ct)  77  (1  -  ^)  , 

V 

where  a  is  an  integer,  the  a^    are  the  zeros  of  the  left-hand  sides  in  (13)  and 
(lU),  and  where  co    is  considered  as  the  independent  variable.     The  fonctions 
a(t)   and  a'(t)  are  independent  of  co  and  will  be  discussed  presently.     For  any 
given  real  t,  a  is  either  one  or  zero  ,  depeinding  on  whether  co  •  0  is  or  is 
not  a  zero  of  the  left-hand  sides  in  (13),   (Hi). 

To  show  that  a(t)   -  a'(t),  we  may  let  a>->ioo.     In  this  case  we  have 
from  Section  IIIj 

(15)  lira      C^  -  1, 

co->ico 

and  from  the  argument  used  in  deriving  statement  3  we  can  show  that 

(16)  lira      e^"'^    y  (n/2,<o,t)  -  1/2   . 
co->ioo 

Therefore,  the  quotients  in  (12)   tend  to  unity  if  a5->ioo.     This  completes  the 
proof  of  equation  (l)   for  real  values  of  t,  and  for  complex  values  of  t     we  ob- 
tain the   same  result  by  analytic  continuation  with  respect  to  t.     The  relation 
(2)   to  (U)  of  theorem  3  can  be  proved  in  the  same  manner. 

We  can  use  Theorem  1  to  compute  the  first  terms  of  the  expansions  of 
C^,  C   ,  S^,  S_^  in  a  series  of  powers  of  t.     If  we  define  "W  (x)  by 

(17)  l^(x)-y^(x)   +  2ioi  y2(x), 
we  find  that 

(18)  OJ(x)  =  '^^(x)  +  t:^3^(x)  +  t^-J^g^^)  *   •    •    •     ' 
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where 

<19)  ^  ,(x)  .  e^i'^^ 

and 

(20)  ^  Jx)   -  -  -^  /    sin     2co  (x-5)     cos  2  4  7n„i(«)  d^  • 

o 

The  explicit  formulas  become  rather  complicated  for  n>2  and  for  arbitrary- 
values  of  X,  but  they  become  simpler  if  x  is  a  miiltiple  of  n/2.     The  results 
derived  from  (l)   to  (7)   and  (17)   to  (19)  are; 

(21)  c  -1.  -fe^  ^^ — (-r-  -"^«'^)*  •  •  ' 

hiJ- -\       oj(Uu  -1)     ^  Uci  -  1 

2 

(22)  S^  -  1  +     ^  (  -^    -    -~-     ctg  na>  )  +   .   .   . 

CO  CO  -1  Uoi    -1 

2 

(23)  C^  -  1  -    ^  f \-      +     — ^  ctg  no>)  +   .   .    . 

co(co  -1)  140)  -1 

2 

(2U)     S_  .  1 1^    *     -^  (  -^ ♦     ^  tan    Tico)  *   .   .    . 

Uco  -1  iico  -1       Uo)   -  1 

The  series  converge  for  all  values  of  t.     It  can  be  shown  that  the 
coefficient  of  t^  in  the  expansion  of  any  one  of  the  expressions  (l)  to  (U)  is  of 
the  type 

R-(o>)  cos  con  +  Rp(co)   sin  con  , 

where  R-,   and  Rp  are  rational  fxmctions  of  co  which  will  have  no  poles  except  for 

CO  »  0,   ±1,...  ,   ±  ('JU  1)   and  (0  -  +  ^  ,  i  I  ,    .   .   .   ,  i  (-n-  i)  . 

Also,   for   lco|->oo   , 

lim     CO     R(co) 

C0->00 

exists,   where  R  denotes  either  R,   or  E_  and  where  k  =  n  +  1  for   the  expression 
in  (3) J   k  =  n  for  the  expressions  in  (l)  and  (U),  and  k  =  n-1  for  the  expression 
in  (2). 
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V.  Fourier  Transfomis  with  Respect  to  the  Parameter  ti>. 

In  this  section  the  notations  of  the  previous  sections  will  be  used» 
We  may  summarize  the  results  as  follows: 

Theorem  U.     Let     y(xja3,t)  be  a  solution  of  Mathieu's  equation  with 
the  initial  VcLlues 

(1)  y(0;a),t)   -  a,  y'(0}w,t)   =  b. 

Then  we  have  for  real  values  of  x  and  od: 

+  X 

(2)  y(x;w,t)  -  a  cos  2cox  +   J     e^^'^  H(x-,^,t)  d^f  , 

-X 

where  for  all  t,  H(x;j^,t)   can  be  expanded  in  a  series  of  powers  of  t.   The  first 
terms  of  the  expansion  are 

(3)  H(x;9f,t)  "     2  +  ^r*2a  sinx  cos0  +  2b cos  x  (cosx  -  cos(2fJ 

.t2[.  ..]. 
In  particular,  if  a  ■  0,  b  «  1,  we  find  that 

(U)      y2(x;o.,t)  -  /  e^^"^  G(x,^,t)  d^), 

-X 

where  G  can  be  expanded  in  a  power  series  in  t  which  is  everywhere  convergent. 
As  a  function  of  the  real  variables  ^  and  x,  G  satisfies  the  partial  differen- 
tial equation 

2     2 

(5)  ^  -  -^  +  Stcos  2x  G  -  0. 


dx     3^ 


3G 


(6)  G(x;ix,t)  »  X  ,  QtJ^vi'it.   x,t)  -It  sinx  cosx  . 

The  proof  of  Theorem  U  can  be  derived  from  Theorem  X  ,  p. 13  in  the 
book  by  Paley  and  Wiener  [^8 J  .  According  to  this  theorem,  the  following  two 
classes  of  functions  are  identical: 

(l)  The  class  of  all  entire  functions  F(z)  satisfying 

(7)  |F(z)|  -  o(e^'^l)  ; 
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(II)  the  class  of  all  entire  functions  of  the  form 
+A 


(8)  F(z)   =/      f(u)   e^^V, 


-A 
where   f(u)   belongs   to  Lp  over   (-A,A), 

It  may  suffice  to  prove   (U),   (5)  and  (6).     Putting 

fn\  t       \         sin  2cox 

(9)  u^(x,co)  -     —^ 

and 

)       ^ 

(10)  u^(x,co)  -  -     ^  y    sin     2a)(x-f)     cos24  u^__j^(4)  d?  , 

o 

we  find 

00 

(n)  72^*'"**^  "     21    ^^n^"  » 

n«0 

where  the  series  on  the  right-hand  side  converges  for  all  values  of  t.     We 

derive  from  (10)   that  for  »->  loo  , 


-n-1 
|u^\.x,co;  I    =      u\.o- 

Therefore 


(12)  luJx,co)|   =     0(0.-"--^)    . 


(13)         g^(x,9^)  .     i/°e-2i'^     u^(x,co) 


exists  and  is  (n-l)  tunes  differentiable  with  respect  to  <)   with  a  continuous 
(n-l)-st  derivative.  Also,  we  find  from  (9)  and  (10)  by  induction  that,  as 
a  function  of  the  complex  variable  co^ 

(lU)  K(x,a.)|  =  o(e2|*"l  )  . 

Therefore,   Paley  and  Viiener's  theorem  proves  that  g  (x,(2()  vanishes  for  real  x, 
9f  if    |x|>  |{2(|,   and  then  we  find  (13)  by  Fourier's  inversion  formxila 

^     2ica5? 


(15)  ■u^(x,a.)  -    /    h"^'^  gjx,i2)  d0  . 
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From  the  differentiability  of  g  {x,(j()   with  respect  to  g   it  also  follows  that  for 
n  >  2 

^^n     ^  ^n 

vanish  for  ^  =  i  x.     By  substituting  the  expressions  in  (15)   for  u     and  u     , 
into  equation  (10)  we   find  after  an  application  of  Fourier's  inversion  formula: 

(17)  g„(x,?)  -  -  ^  /V-i^-io.  /lisM-S)  cos  2K  d«    /V»»''g„.,(4,r)dr 

-00  -4 

Exchanging  the  order  of  integrations  and  carrying  out  the  integration  with  res- 
pect to  CO  we  obtain  the  following  result:  Let 

(18)  L(x,4,^  )  =  / "^  e^^"^^  -^^  sin  2a)(x-4)  f  ; 

-00 

then 

X     E 

(19)  g^U,9)   -  -  ^  /  d4  /d>^  L(x,4,?,r)  gn-1^^'^^  ^°^  ^^   ' 

0    K 

If  the  integrals  from  -co  to   +oo  in  (17)  and  (19)   are  defined  as 

A 
lim     y  , 
A->co   -A 

the  only  possible  values  for  L  are  0,1  s-n  ,  *n  .     Since  g     «  •^  ,  it  follQvs 

by  induction  from  (19)   that  for  -x  '^  91  ^  x 

lg^(x,9)  <  2^"-^   |x|2"  /  nl 

Therefore, 

'  00 

G(x,9i,t)   -    ^    g„(x,^)t" 

n=0 

exists  for  all  values  of  t.     That  G  satisfies  (5)>subject  to   the  initial  con- 
diticriE  (6),  follows  from  integration  by  parts  of 


/  e^i'^ 
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i^  _  a_G  ^  Qtcos  2xG  \d(^ 


) 


and  from  equation  (lO)  which  implies  that 


,2 

n      2 
K    +  l4^  u  ■  -8t  (cos  2x)  u  , 

dx 


VI.  Generalizations ♦ 

Most  of  the  results  of  the  previous  sections  can  be  generalized 
very  easily.  Let 

(1)  y*  +  Ilw^  +  8f(x)  ]  y  -  0 

be  a  differential  equation  of  Hill's  type.     In  particular,  we  shall  assume  that 

00 

(2)  f(x)  '    Yl    ^     cos2nx 

n=l    ^ 

is  an  even  function  of  period  n   which  we  shall  assume  to  be  analytic  everywhere. 
Then  we  can  define  the  analogues  of  the  infinite  determinants  A  (co.t),  C  ,  S  , 

0+    + 

C  ,  S  and  vb   can  prove  all  the  inequalities  in  Section  III  if  we  replace  t  by 

]^  |t  I  .  If  we  define  the  standard  solutions  y,  and  y^  by  the  same  initial 

conditions  as  before,  we  find  exactly  the  same  conditions  under  which  (1)  has  a 
periodic  even  or  odd  solution  of  period  n  or  2Tijand  we  can  prove  (lV.5),  (IV. 6) 
and(IV.7).  We  can  also  prove  the  results  of  section  V  after  an  appropriate  modi- 
ficAion  of  (V.5)  and  (V.6). 

But  a  difficulty  arises  if  vre  wish  to  generalize  the  equations  (1) 
to  (U)  of  Section  IV.  To  do  this  we  would  have  to  prove  that  the  right-hand 
sides  in  (IV.I)  to  (IV .U),  or  rather  their  generalizations  in  the  present 
case,  have  simple  zeros  with  respect  to  u>  if  the  t  are  real.  We  would  need 
either  the  analogues  of  the  product  formulas  (II.2U)  and  (11,25)  for  the  corres- 
ponding infinite  determinants  or  a  substitute  for  them.  So  far  this  difficulty 
has  not  been  overcome . 
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VII.  Application  to  the  Construction  of  Transparent  Media 


Let  f(x)  be  a  continuous  function  of  x  such  that 
if  X  <  0  or  X  >  Z  ,  and 


(1)  f(x)  =  k^ 


(2)  f(x)  =  y  c  cos  2mx 

*—   n 
m=o  - 

for  0  =  X  ^  X  where  >C  =  nn  ,  n  =  1,2,3,  ....  We  may  consider  ^f(x)^  as  the 
index  of  refraction  of  a  mediuin  which  is  homogeneous  except  for  a  layer  of 
thickness  -C  .  A  wave  traveling  in  this  mediuin  will  be  defined  by  a  function 
y(x)  which  satisfies 

(3)  y*  +  f(x)y-  0 

and  is  continuous  and  has  a  continuous  first  derivative , 

We  shall  call  the  medium  completely  transparent  if  for  any  real  value  of 
k  there  exists  a  solution  y(x)  as  defined  above  which  satisfies 

(U)  y(x)   =  e^^  for  x  <  0,  y(x)   =  e^''^''"'"^   forx>i  , 

•5t- 

where  a  is  real.  . 

Let  y, (x)  and  y2(x)  be  the  standard  solutions  of  (3)  in  0$x  =  /  , 

which  are  defined  by  y^(0)   -  y2(0)   -  1,  y2(0)  -  7-^(0)   -  0    . 

Then,  for  0  =  x  <  K 

(5)  y(x)  =  y-^U)  +  iky2(x) 

because  y(x),  y  (x)  are  continuous  for  x  «  0.  For  x  »  "we  find  from  (U)  : 

(6)  ly(/)|  -  1,        y'(i)  -  iky(/), 
and  therefore 

(7)  yld)   +  k^  y|(/)  =1,   y^ii)  '  ikyJ^Jt)  ,  y[d)  -y^ii). 


■«•  Dr.  J.  Shmoys  gave  the  following  physical  interpretation  to  this  definition: 
If  the  layer  is  submerged  in  a  mediuin  of  the  constant  index  of  refraction 
k^,  a  plane  wave  passes  through  the  layer  without  reflection,  independent 

of  the  value  of  k^. 
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If  (7)  is  true  for  all  k  then 

(7')  y-L(^)  «=  1  1,      j^il)  '0,  y^d)  "0,   ygC/)  -  1   1. 

This  shows  that  the  following  is  a  necessary  and  sufficient  condition  for  a 
transparent  layer:   Tl-ie  differential  equation  (3),  in  which  f(x)  is  now  given 
by  (2)  for  all  values  of  x,  must  have  two  linearly  independent  solutions  having 
periods  of  either  t  or  2 /6  . 

The  remarks  in  Section  II  show  that  hathieu's  equation  can  indeed  be 
used  for  the  construction  of  transparent  layers  end   that  this  is  also  true  in 
general  for  any  equation  of  Hill's  type  (3)  if  we  choose  /=  3n,  Un,  ^n,    .... 
We  merely  have  to  make  sure  that  the  characteristic  exponent  a   of  Section  II 
is  a  rational  number  between  0  and  ^  .  This  has  been  pointed  out  already  by 
Kramers  TUl .  If  we  take  k  ■  f(o)  »  f(-c  ),  we  obtain  an  index  of  refraction 
which  is  continuous  and  has  a  continuous  first  derivative  everywhere.  But  in 
the  cases  obtained  from  Kathieu's  equation,  the  inhomogeneous  part  of  the  layer 
consists  of  at  least  three  equal  sublayers.  Cases  where  -0  ■  n,  that  is  where  the 
period  of  f(x)  itself  is  also  a  period  of  all  solutions  of  (3)^  have  been  found 
by  Meissner  f6]  ,  who  used  a  discontinuous  piecewise  constant  f (x),  and  by 
Klotter  and  Kotowski  [31  ,  whose  method  consists  in  a  numerical  analysis  of 
the  regions  of  stability  for  the  solutions  of 

(8)  y  +  (X  +  Y^  cosx  +  Yp  cos2x  )  y  »  0. 

In  [33 two  periodic  solutions  of  period  n  are  then  found  for  y-i  =  1>  Yp  =  V^ 

and  \**'1.17.  Finally,  Bouwkamp  [I'}   has  shown  that  the  case  where  f(x)  repre- 
sents a  periodically  repeated  triangle  can  be  discussed  completely. 

f  1  '/* 

If  one  wishes  to  interpret  |f(x) }    as  the  index  of  refraction  of 

a  transparent  layer,  one  may  i-eqiire  an  f(x)  which  is  positive  everywhere.  It 

is  easily  verified  that  the  minima  of 

(9)  f(x)  e  X  +  Yi  cosx  +  Y2COE2  X 
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are  at  the  points  at  which  either  sinx  «  0  or  cosx  «  ~Y2^/(Uy2)    •     ^'^'^  '^^^ 
values  of  y-,>Y2*^  ^s®^  i^f^D*  "this  gives,   at  the  minimum,  i.e.,  x  =  5/(6'n), 

(10)  f  (|iT  )  -  ^  -  ^  Yi  -  I  Y2     >    0 

and  at  X  "  Ti 

(11)  tU)  =  X  -  Yi  >  0- 

Therefore,  the  result  obtained  in  ["3  3  provides  such  an  example. 

A  last  remai'k  to  be  made  in  connection  vith  transpaj-ent  layers  is 
the  following  one:  I^et 

(12)  k^  =  f(o)  =  f(n)  «  f(£) 

and  consider  the  phase  9>(x)  of  the  incoming  wave  exp(ikx)  for  (x  <  O)  in  the 
layer*  We  can  define  (jSix)   by 


o  y^^C^)  +  k  j^iK) 
since  it  can  be  shown  that 

1 
(lU)  y^^Cx)  +  ikygCx)  -(y^(x)  +  y|(x)|  ^  exp(i9f(x)  )  . 

(This  is  an  elementary  consequence  of  the  fact  that  y,   and  y^  satisfy  a  homo- 
geneous linear  differential  equation  of  second  order  which  does  not  have  a  term 
involving  y   .)     The  increase  in  phase  of  the  wave  at  the  end  of  the  layei'  will  be 

(15)     g(   )  -  !?(n)  «  2nTi  ,  (n  «  1,2,3,    ...   ) 

since  the  left-hand  side  in  (lU)  must  be  unity  for  x  »  n  according  to  (7').     It 
is  easily  seen  that  the  integer  n  does  not  deyiend  on  k,  because  the  indefinite 
integral  on  the  right-han^   side  in  (13)   can  be  written  as 

ky2(x) 

Therefore,  ^(n)  depends  only  on  the  number  of  zeros  of  y2(x)  in  the  interval. 
0  «  X  ■$  Ti.  Without  the  layer  the  phase  of  the  incondng  wave  at  x  »  n  would 
have  been  kn  .  If  it  would  be  possible  to  choose  k  «{  f(o,.l]    in  such  a  way 
that  k  ■  2n,  where  n  is  the  same  as  in  (15),  we  would  have  a  layer  which  would 
not  even  cause  a  loss  of  phase.  So  far  a  construction  of  such  a  layer  has  not 
been  carried  out,  but  uhere  does  not  appear  to  be  any  reason  why  it  cou!l-d  not 
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exist.     Of  course,   the   tJri%d.al  case  where  f(x)  ;  k     shoiild  not  be  considered. 
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